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Abstract. The concordance Cold Dark Matter model for the formation of structure in the 
Universe, while remarkably successful at describing observations on large scales, has a number 
of problems on galaxy scales. The Milky Way and its satellite system provide a key laboratory 
for exploring dark matter (DM) in this regime, but some of the most definitive tests of local DM 
await microarcsecond astrometry, such as will be delivered by the Space Interferometry Mission 
(SIM Planetquest). I discuss several tests of Galactic DM enabled by future microarcsecond 
astrometry. 
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1. Introduction 

Since the seminal study of Searle & Zinn (1978) the notion of accretion of "subgalactic 
units" , including "late infall" , has been a central concept of stellar populations studies. 
N-body simulations of the formation of structure in the Universe in the presence of 
dark matter (and dark energy) also show galaxies (and all large structures) building up 
hierarchically. However, while the active merging history on all size scales demonstrated 
by high resolution. Cold Dark Matter (CDM) numerical simulations have had remarkable 
success in matching the observed properties of the largest structures in the Universe 
(like galaxy clusters), they are a challenge to reconcile with the observed properties 
of structures on galactic scales. The Milky Way (MW) and its satellite system are a 
particularly important laboratory for testing specific predictions of the CDM models. 
The era of microarcsecond astrometry such as will be delivered by SIM Planetquest, will 
enable a number of definitive dynamical tests of local CDM by way of determining the 
distribution of Galactic dark matter (DM). We discuss several of these tests here. 

2. Measuring 6lsr with the Sagittarius Stream 

An understanding of the distribution and amount of DM in the Galaxy — indeed, 
virtually every problem in Galactic dynamics — depends on establishing the benchmark 
parameters _Roj the solar Galactocentric distance, and OlsR: the Local Standard of Rest 
(LSR) velocity. A 3% error in both i?o and 8lsr leads to a 5% error in the determination 
of the Galactic mass scale using, e.g., traditional Jean's equation methods. 

Despite decades of effort, the rate of Galactic rotation at the Sun's position remains 
uncertain, with measurements varying by 20% or more. Hipparcos proper motions have 
been used to determine that 8lsr = (217.5 ± 7.0)(i?o/8) km s^^ using disk Cepheid 
variables (Feast k Whitelock 1997) and (240.5 ± 8.3)(i?o/8) km s^^ from OB stars 
(Uemura et al. 2000), while Hipparcos proper motions of open clusters have yielded 
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207(i?o/8) km (Dias & Lepine 2005) [though see the alternative analysis of these open 
cluster by Princhaboy 2006 and Frinchaboy & Majewski 2006, which yields (221 ± 3) km 
s""'^ for any 7 ^ Rq ^9]. Meanwhile, Hubble Space Telescope measurements of the proper 
motions of bulge stars against background galaxies in the field of the globular cluster M4 
yield Olsr = (202.4 ± 20.8) (J?o/8) km s'^ (Kalirai et al. 2004) and (220.8 ± 13.6) (i?o/8) 
km s"-'^ (Bcdin ct al. 2003). Radio measures of the proper motion of Sgr A* (Rcid & 
Brunthaler 2004) yield a transverse motion of (235.6 ± 1.2)(i?Q/8) km s~^, but, when 
corrected for the solar peculiar motion, yield ©lsr = 220(iZo/8) km s~^ (M. Reid, this 
proceedings). Of course, as may be seen, most of these measures depend on an accurate 
measure of Rq as well as knowledge of the solar motion (which is apparently now more 
controversial than previously suspected; cf. the significantly smaller solar motion derived 
by Dehnen & Binncy 1998 compared to that from M. Reid, this proceedings). Moreover, 
considerations of non-axisymmetry of the disk yield corrections to the measurements 
that suggest Olsr may be as low as 184 ± 8 km (Oiling & Merrifield 1998) or lower 
(Kuijken & Tremaine 1994). 

Based on recent estimates of the luminosity of the MW (Flynn et al. 2006), any Olsr 
exceeding 220 km s~^ puts our Galaxy more than la away from the TuUy-Fisher relation. 
Clearly, additional independent methods to ascertain ©lsr are valuable for establishing 
with certainty the MW mass scale and whether our home galaxy is unusual. 

An ideal method for ascertaining ©lsr would be to measure the solar motion with 
respect to a nearby reference known to be at rest in the Galactic reference frame (at least 
in the disk rotation direction), since this would not only overcome traditional difiiculties 
with working in the highly dust-obscured and crowded Galactic center, but be indepen- 
dent of Rq and any assumptions that the reference lies in the center of the MW potential. 
The debris stream from the tidally disrupting Sagittarius (Sgr) dSph galaxy provides al- 
most nearly this ideal situation. This dSph and its extended tidal debris system orbit the 
MW in almost a polar orbit, with a Galactic plane line-of- nodes almost exactly along the 
Galactic X-axis. The remarkable coincidence that the Sun presently lies within a kilo- 
parsec of the Sgr debris plane, which has a pole at {Ip, bp) = (272, —12)° (Majewski et al. 
2003), means that the motions of Sgr stars within this plane are almost entirely contained 
in their Galactic U and W velocity components, whereas their V motions almost entirely 
refiect solar motion. Thus, in principle, the solar motion can be derived directly by the 
m cos(6) dimension of the Sgr stream proper motion — though refined results come from 
comparing to Sgr models that can account for deviations from the ideal case (Majewski 
ct al. 2006). An additional advantage is that stars in the Sgr stream, particularly its 
substantial M giant population (Majewski et al. 2003), are ideally placed for uncrowded 
field astrometry at high MW latitudes, and at relatively bright magnitudes. 

This combination of factors means this experiment requires only the most modest pre- 
cisions from SIM, is well within the capabilities of Gaia, and may even be within the grasp 
of future high quality, ground-based astrometric studies. The latter was demonstrated 
by ground-based astrometry obtained by Casetti-Dinescu et al. (2006) in two Kapteyn 
Selected Area fields positioned on the Sgr trailing arm, where dozens of 1-3 mas year~^ 
proper motions for V = 19 — 21 main sequence stars in the Sgr stream were averaged to 
derive a mean motion; these data showed that with more stars and more fields at this 
level of astrometric precision one might hope to distinguish MW models at the level of 
10s of km s~^. To achieve ~ 1 km s~^ precision in the solar motion with this method 
will require ~ 0.01 mas year~^ mean proper motions for sections of the Sgr trailing arm 
— a result that will be trivial for microarcsecond level, space astrometry. 
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3. Probing the Galactic Potential with Tidal Streams 

It is now well known that the Galactic halo is inhomogeneous and coursed by streams 
of debris tidally pulled out of accreted satellites. Standard methods of measuring the 
MW's gravitational potential using a tracer population whose orbits are assumed to be 
random and well-mixed are systematically biased under these circumstances (Yencho et 
al. 2006). However, these streams themselves provide uniquely sensitive probes of the 
MW potential. If we could measure the distances, angular positions, radial velocities and 
proper motions of debris stars, we could integrate their orbits backwards in some assumed 
Galactic potential. Only in the correct potential will the path of the stream stars ever 
coincide in time, position and velocity with that of the parent satellite (see Fig. 1). 

Microarcsecond proper motions and parallaxes combined with ground-based radial 
velocities provide everything needed to undertake the experiment. If we find a coherent 
stream without an associated satellite, the same techniques apply, but with the parent 
satellite's position and velocity as additional free parameters. The usefulness of Galactic 
tidal streams for probing the MW potential has long been recognized, and the promise of 
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Figure 1. A demonstration of the use of tidal streams to measure tlie sliape and strength of 
the Galactic potential if precise 6-D information is available for stream stars. A Sgr-like tidal 
stream was created by the disruption of a dwaxf satellite in a rigid Galactic potential through 
a semi-analytical N-body simulation (e.g., as in Law et al. 2005; panel a). With complete 6-D 

phase space information in hand, guesses may be made on the strength and shape of the Galactic 
potential, and the orbits of the individual stars in the tidal streams run backwards under these 
assumed potentials. Panels (b) and (c) demonstrate what happens when the strength of the 
Galactic potential is underestimated by varying degrees: When the orbits are run backwards, 
the tidal stream stars orbit at too large a radius and do not converge to a common phase space 
position. In panels (e) and (f ) the strength of the Galactic potential has been overestimated, and 
when the clock is run backwards the tidal stream stars assume orbits that are too small and once 
again do not converge on a common phase space position. In panel (d) a Galactic potential of 
the correct strength was guessed, and, when the stream star orbits are run backwards, the tidal 
stream stars collect back into the core of the parent satellite (shown as the dot in all panels). 
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astrometric space missions to provide the last two dimensions of phase space information 
for stream stars eagerly anticipated (e.g., Johnston et al. 1999, Penarrubia et al. 2006). 

Applying this method to simulated data observed with the /ias yr"'^ precision proper- 
motions possible with SIM and km s~^ radial velocities suggests that 1% accuracies on 
Galactic parameters (such as the flattening of the potential and circular speed at the 
Solar Circle) can be achieved with tidal tail samples as small as 100 stars (Johnston 
et al. 1999, Majewski et al. 2006). Dynamical friction is not an important additional 
consideration if the change in the energy of the satellite's orbit in Norh orbits is less than 
the range in the energies of debris particles. For A'orb = 3, this condition is met for all 
satellites except the Magellanic Clouds and Sgr. Evolution of the Galactic potential does 
not affect the current positions of tidal debris, which respond adiabatically to changes 
in the potential and therefore yield direct information on the present Galactic mass 
distribution independent of how it grew (Pefiarrubia et al. 2006). 

It is important to sample tidal tails to points where stars were torn from the satellite 
at least one radial orbit ago and have thereby experienced the full range of Galactic 
potential along the orbit (Johnston 2001). The stars also need to have proper motions 
measured sufficiently accurately that the difference between their own and their parent 
satellite's orbits are detectable — this translates to requiring proper motions of order 
100 /Ltas yr~^ for the Sgr stream but 10s of /zas yr~^ for satellites farther away. Ideally, we 
would like to probe stars from different tidal tails at a variety of Galactocentric distances 
and orientations with respect to the Galactic disk to derive the shape and density profile 
of the DM. In the last few years evidence for dSph galaxy tidal tails has been discovered 
around systems at large Galactic radii (Munoz et al. 2006a, Sohn et al. 2007) and, with 
SIM, extended versions of such distant streams can be used to trace the Galactic mass 
distribution as far out as the virial radius with an unprecedented level of detail and 
accuracy. 

A key problem that is currently debated but that will easily be resolved with precision 
astrometry of tidal tails is the shape of the Galactic halo. CDM predicts that host DM 
halos should be oblate (minor to major axis ratio, q ~ 0.7) with flattening increasing 
with radius, whereas warm and self-interacting DM models tend to find rounder halos 
(Bullock et al. 2002). Studies of the halos of external galaxies — whose shapes can be 
measured via flaring of HI disks, the shapes of diffuse x-ray isophotes (for ellipticals), 
or with polar rings — typically find rather oblate results (see summaries in Combes 
2002, Merrifield 2005). Meanwhile a controversy has emerged regarding the true shape 
of the halo as measured by the Milky Way's own polar ring system — the tidal stream 
of the Sgr dwarf spheroidal galaxy. Depending on how the extant data are analyzed and 
modeled the debris of the Sgr dwarf has been argued to be consistent with a rather 
prolate {q ~ 1.4; Helmi 2004), near spherical (Ibata et al. 2001, Fellhauer et al. 2006), 
slightly oblate (g ~ 0.93; Johnston et al. 2005, Law et al. 2005) or oblate {q = 0.85; 
Martinez-Delgado et al. 2004) Galactic halo. The definitive measurement of the shape 
of the Galactic halo within the Sgr orbit will come from SIM proper motions of stars 
in its tails, and the variation of the halo shape with radius will come by repeating the 
technique with other tidal streams at other distances. Interestingly, recent measures of 
the mass profile of the MW using stellar tracers suggest that it gets rounder with radius 
(Kuijken 2003), in contradiction to the predictions of CDM. 

The possible existence of a significant fraction of the halo in the form of dark satellites 
has been discussed in recent years (Moore et al. 1999, Klypin et al. 1999). These putative 
dark subhalos could scatter stars in tidal tails, possibly compromising their use as large- 
scale potential probes, but astrometric measurements of stars in these tails could, on the 
other hand, be used to assess the importance of dark substructure (Ibata et al. 2002, 
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Johnston et al. 2002). Early tests of such scattering using only radial velocities of the 
Sgr stream suggest a MW halo smoother than predicted (Majewski et al. 2004), but this 
represents debris from a satellite with an already sizable intrinsic velocity dispersion. 
Because scattering from subhalos should be most obvious on the narrowest, coldest tails 
(such as those seen from globular clusters, like Palomar 5 - e.g., Grillmair & Dionatos 
2006) these could be used to probe the DM substructures, whereas the stars in tails 
of satellite galaxies such as Sgr, with larger dispersions initially and so less obviously 
affected, can still be used as global probes of the Galactic potential. 

Gaia will allow these tests to be attempted for nearby streams, whereas with SIM, 
it will be possible for the first time to probe the full three-dimensional shape, density 
profile, extent of and substructure within our closest large DM halo. 

4. Probing the Galactic Potential with Hypervelocity Stars 

A complementary method for sensing the shape of the Galactic potential can come from 
SIM observations of hypervelocity stars (Gnedin et al. 2005). Hills (1988) postulated that 
such stars would be ejected at speeds exceeding 1000 km s~^ after the disruption of a close 
binary star system deep in the potential well of a massive black hole, but HVSs can also 
be produced by the interaction of a single star with a binary black hole (Yu & Tremaine 
2003). Recently Brown et al. (2006) report on five stars with Galactocentric velocities 
of 550 to 720 km s""'^, and argued persuasively that these must be "imboimd stars with 
an extreme velocity that can be explained only by dynamical ejection associated with a 
massive black hole" . After the success of these initial surveys, it is likely that many more 
HVSs will bo discovered in the next few years. 

If these stars indeed come from the Galactic center, the orbits are tightly constrained 
by knowing their point of origin. In this case, as Gnedin et al. (2005) demonstrate, 
the non-spherical shape of the Galactic potential due in part to the flattened disk 
and in part to the triaxial dark halo will inchicc non-radial inflections (which will be 
primarily in transverse direction at large radii) in the velocities of the HVSs of order 5-10 
km s^^, which corresponds to 10 100 /ias yr^^. Each HVS thus provides an independent 
constraint on the potential, as well as on the solar circular speed and distance from the 
Galactic center. The magnitudes of the known HVSs range from 16 to 20, so their proper 
motions will be measurable by SIM with an accuracy of a few ^as yr~^, which should 
define the orientation of their velocity vectors to better than 1%. With a precision of 20 
fias yr~^ the orientation of the triaxial halo could be well-constrained and at 10 /uas yr~-^ 
the axial ratios will be well- constrained (Gnedin et al. 2005). 

5. Dark Matter within Dwarf Galaxies 

Dwarf galaxies, and particularly dSph galaxies, are the most DM-dominated systems 
known to exist. Relatively nearby Galactic dSph satellites provide the opportunity to 
study the structure of DM halos on the smallest scales, and, with microarcsecond as- 
trometry, make possible a new approach to determining the physical nature of DM (Stri- 
gari et al. 2007, 2008). CDM particles have negligible velocity dispersion and very large 
central phase-space density, which results in cuspy density profiles over observable scales 
(Navarro et al. 1997, Moore et al. 1998), whereas, in contrast. Warm Dark Matter (WDM) 
has smaller central phase-space density so that constant central cores develop in the den- 
sity profiles. Because of their small size, if dSph cores are a result of DM physics then 
the cores occupy a large fraction of the virial radii, which makes these particular cores 
more observationally accessible than those in any other galaxy type. Using dSph central 
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velocity dispersions, earlier constraints on dSph cores have excluded extremely warm 
DM, such as standard massive neutrinos (Lin & Faber 1983, Gerhard & Spergel 1992). 
More recent studies of the Fornax dSph provide strong constraints on the properties of 
sterile neutrino DM (Goerdt et al. 2006, Strigari et al. 2006). 

The past decade has seen substantial progress in measuring radial velocities for large 
numbers of stars in nearby dSph galaxies, and the projected radial velocity dispersion 
profiles are generally found to be roughly flat as far out as they can be followed (Muhoz 
et al. 2005, 2006a, Walker et al. 2006a,b, 2007, Sohn et al. 2007, Koch et al. 2007a,b, 
Mateo et al. 2007). By combining such radial velocity profiles with the surface density 
distributions of dSph stars (which are well-fitted by King profiles, modulo slight varia- 
tions, especially at large radii), it is typical to derive DM density profiles by assuming 
dynamical equilibrium and solving the Jeans equation (e.g, Richstone & Tremaine 1986). 
The results of such equilibrium analyses typically imply at least an order of magnitude 
more dSph mass in dark matter than in stars as well as mass-luminosity ratios that 
increase with radius — in some cases quite substantially (e.g., Kleyna et al. 2002) — 
though at large radii tidal effects probably complicate this picture (Kuhn 1993, Kroupa 
1997, Munoz et al. 2006a, 2007, Sohn et al. 2007, Mateo et al. 2007), and in some cases 
mass-follows-light models incorporating tidal disruption describe the observations quite 
well (Sohn et al. 2007, Mufioz et al. 2007). Knowing whether mass follows light in dSphs 
or if the luminous components lie within large extended halos is critical to establishing 
the regulatory mechanisms that inhibit the formation of galaxies in all subhalos. 

Unfortunately, as shown by Strigari et al. (2007), equilibrium model solutions to the 
Jeans equation are degenerate and satisfied by DM density profiles with either cores or 
cusps. This is because there is a strong degeneracy between the inner slope of the DM 
density profile and the velocity anisotropy, /3, of the stellar orbits, which leads to a strong 
dependency of the derived dSph masses on f3. Radial velocities alone cannot break this 
degeneracy (Fig. 2) even if the radial velocity samplers htq increased to several thousand 
stars (Strigari et al. 2007). The problem is further compounded if we add triaxiality. 
Galactic substructure, and dSph orbital shapes to the allowable range of parameters. 
The only way to break th(^ mass-anisotropy degeneracy is to mcasiire more phase space 
coordinates per star, and in particular to acquire transverse velocities for the stars, 
because the Jeans equation solved in the transverse dimension probes the anisotropy 
differently than in the radial velocity dimension. Thus, combining high precision proper 
motions with even the present samples of radial velocities holds the prospect to break 
the anisotropy-inner slope degeneracy (Fig. 2). 

The most promising dSphs to obtain proper motions for will be the nearby (60-90 
kpc distant) systems Sculptor, Draco, Ursa Minor, Sextans and Bootes, which include 
the most DM-dominated systems known (Mateo 1998, Mufioz et al. 2006b, Martin et al. 
2007) as well as a system with a more modest M/L (Sculptor). To sample the velocity 
dispersions properly will require proper motions of > 100 stars per galaxy with accuracies 
of 7 km s^^ or better (less than 15 /las yr~^). Strigari et al. (2007) show that with about 
200 radial velocities and 200 transverse velocities of this precision it will be possible to 
reduce the error on the log-slope of the dark matter density profile to about 0.1 — which 
is an order of magnitude smaller than the errors attainable from a sample of 1000 radial 
velocities alone, and sensitive enough to rule out nearly all WDM models (Fig. 2). Figure 
3 shows that exploring even these nearby dSph systems requires precision proper motions 
of stars to y > 19, a task well beyond the capabilities of GAIA, but well-matched to the 
projected performance of SIM, though requiring a Key Project level of observing time 
(e.g., 100 days of SIM observing for 200 Draco stars). Note that while the Sgr dSph is 
several times closer than these other dSphs, it is obviously a system in dynamical non- 
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Figure 2. Left panels: A demonstration of the ability to recover information on the nature of 
DM using observations of dSph stars, from analytical modehng by Strigari et al. (2007) . ElUpses 
indicate the 68% and 95% confidence regions for the errors in the measured dark halo density 
profile slope (measured at twice the King core radius) and velocity anisotropy parameter /3 in the 
case where only radial velocities are available for 1000 stars in a particular dSph (top panels). 
A significant improvement is derived from the addition of 200 SIM proper motions providing 5 
km s~^ precision transverse velocities (bottom panels). The left (right) panels correspond to a 
cusp (core) halo model for dSphs and the small x's indicate the fiducial, input model values. 

Figure 3. Right panel: Potential SIM PlanetQuest exploration of the Draco dSph would need 
to probe \,oV ^ M = 19 to derive a sample of 200 red giants as seen by its Washington Af-band 
luminosity function {thick line and left axis). The thin lines represent the number of days {right 
axis) necessary to observe all the stars to a given magnitude limit with SIM and for a given 
transverse velocity uncertainty: 3 km s~^ {thin solid line), 5 km s~^ {dotted line), 7 km s~^ 
{dashed line) and 10 km s~^ {dotted-dashed line). Prom Strigari et al. (2008). 



equilibrium; exploring this system will be of interest as a case study for establishing the 
internal dynamical effects of tidal interaction on dSphs. 
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